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We treat the non-relativistic Compton scattering process in which an incoming photon scatters 
from an N-electron many-body state to yield an outgoing photon and a recoil electron, without 
invoking the commonly used frameworks of either the impulse approximation (lA) or the independent 
particle model (IPM). An expression for the associated triple differential scattering cross section is 
obtained in terms of Dyson orbitals, which give the overlap amplitudes between the A''-electron 
initial state and the (A'' — 1) electron singly ionized quantum states of the target. We show how 
in the high energy transfer regime, one can recover from our general formalism the standard lA 
based formula for the cross section which involves the ground state electron momentum density 
(EMD) of the initial state. Our formalism will permit the analysis and interpretation of electronic 
transitions in correlated electron systems via inelastic x-ray scattering (IXS) spectroscopy beyond 
the constraints of the lA and the IPM. 
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I. INTRODUCTION 

Compton scattering is unique among spectroscopic 
techniques in that it allows direct experimental access 
to the ground state electron momentum density (EMD) 
p(p) of the target many body systemi*^. Recent high 
resolution Compton scattering studies using synchrotron 
light sources have revealed interesting electron correla- 
tion effects in a number of materials^iiASiiifiiSiifliiiii^. 
The experimental work has been concentrated largely 
on the measurement of the double differential scattering 
cross section (for detecting energy transfer and solid an- 
gle of the outgoing photon), which yields the so-called 
Compton profile (CP) related to the two-dimensional 
(2D) integral of the EMD 



J{Pz) = P{p)dpxdpy, 



(1) 



or equivalently a one-dimensional projection of the EMD 
along the direction of the scattering vector pz of the in- 
cident photon. 

Form (QJ which is used in much of the existing analy- 
sis of CP's is obtained within the framework of the im- 
pulse approximation (IA)i2iid. The fundamental scat- 
tering process considered in the lA is the scattering of 
a photon from a collection of free electrons. The I A is 
expected to be valid when the energy transferred in the 
scattering process is much larger than the binding en- 
ergy of the electronic states involved. By its very nature, 
Eq. 1^ lacks a systematic way of taking account of de- 
viations from the lA. With this motivation, our purpose 
in this article is to consider the general scattering event 
in which the incoming photon is scattered from a bound 
many-electron system. We evaluate the resulting par- 
tial triple differential scattering cross section rigorously 
in terms of the so-called Dyson orbitals, which involve 
overlap of the A^-body initial state wave function with 



the (A^ — 1) body wave function of the singly ionized 
final state with an ejected electron. The physically rel- 
evant triple differential scattering cross section is then 
obtained by summing over final (ionic) states and the 
steps necessary to recover the lA are clarified. By going 
beyond the lA, our study provides a systematic scheme 
for understanding electronic structure and correlation ef- 
fects via inelastic x-ray scattering (IXS), away from the 
deeply inelastic regime. 

In this connection, it is important to recognize that 
the standard Compton scattering experiment does not 
involve the measurement of the kinematics of the out- 
going (recoil) electron. This is the reason for the ap- 
pearance of the 2D integral and the concomitant loss of 
information about p(p) in Eq. 1^. As was pointed out 
first by Kaplan and Yudi»l^ on the basis of their theoret- 
ical studiesi^ of Compton scattering on bound electrons 
of light atoms and molecules, the full three-dimensional 
(3D) EMD can be determined, if the characteristics of 
the scattered photon and the ejected electron are mea- 
sured in coincidence The authorsiSiifi have also shown 
that, if the ejected electrons are selected by energy, the 
EMD associated with individual electronic states can in 
principle be obtained. 

Although coincidence experiments were undertaken 
quite earlyiiii^, results for 3D EMD were first reported 
by Bell and collaborator a^'?i^fti^^i^^i^i?i^"' ; see also related 
work of Itoh and collaborator a^'"?'^^ . Since the cross sec- 
tion for an incident photon to scatter into an outgoing 
electron and a photon is measured, such a measurement is 
often referred to as a (7,67) experiment. From a formal 
viewpoint, the (7, 67) experiment provides a measure- 
ment of the triple differential scattering cross section, for 
which we present in this article a rigorous many-body 
expression . 

For interpreting experimental CP's using the I A based 
formula , actual computations in the literature largely 
employ the independent particle model (IPM). The 
many-electron wave function underlying the IPM is built 
from Slater determinants of single-electron orbitals ob- 
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tained usually via unrestricted Hartree-Fock (UHF) ap- 
proach or various versions of the density functional the- 
ory (DFT). With a growing interest in applying IXS for 
investigating electronic transitions in highly correlated 
systems using synchrotron light sources, we should keep 
in mind that one will need to take account of deviations 
not only from the I A but also from the IPM. Even in the 
relatively simple case of Li, substantial deviations in the 
EMD predicted by the local density approximation have 
been implicated in explaining the observed discrepancies 
between the computed and measured CP's^''. 

Concerning other relevant literature related to the is- 
sue of going beyond the lA, several studies have consid- 
ered the accuracy of the lA in describing core Compton 
profile a^^i^^i^°i^^i^^i^^ . A general method for introduc- 
ing final-state-interaction effects has been discussed by 
Sears^i in the context of deep-inelastic neutron scatter- 
ing. This work also discusses the Bjorken-scaling and 
y-scaling properties of the lA which have been particu- 
larly useful in particle physics. Recently, high resolution 
valence Compton profiles (CP's) of Li at a relatively low 
photon energy of 8-9 KeV were considered in Refsi ^^i^^ . 
The observed asymmetries in shape and smearing of the 
Fermi surface features in the CP's were attributed to the 
breakdown of the lA. It is further shown that these dis- 
crepancies in Li can be understood in terms of a finite 
width of the final state spectral function^. To our knowl- 
edge, all previous work concerning the breakdown of the 
lA has been dedicated to understanding the double differ- 
ential scattering cross section. The present study focuses 
on the elementary (7, 67) scattering process and provides 
a clearer picture of the many-body effects and their con- 
nection with the lA and the IPM. 

An outline of this article is as follows. These introduc- 
tory remarks are followed in Section ITlAI with a rigorous 
treatment of the partial triple differential scattering cross 
section in terms of Dyson orbitals. Section lll Bl addresses 
the question of summing over final states to obtain the 
total triple differential cross section and how it reduces 
to the lA result proportional to the EMD. Section III 
makes a few concluding remarks. 

II. GENERAL EXPRESSIONS FOR TRIPLE 
DIFFERENTIAL SCATTERING CROSS SECTION 

A. Partial triple differential cross section 

We consider the non-relativistic elementary scattering 
process in which the incoming photon of energy uji (here 
and throughout this article, natural units, h — c — 1, 
are assumed implicitly) and momentum ki scatters from 
the N-electron many body ground state of the solid with 
energy Eo{N). The final state consists of: an outgoing 
photon with energy uj2 and momentum k2; the {N — 1) 
electron ionized state of the solid characterized by quan- 
tum number n and energy E„(N — 1); a recoil electron 

carrying kinetic energy i?e"'' and momentum p„. We as- 



sume that the momentum transferred to the ionic system 
is q and the associated kinetic energy is neglected given 
the large mass of the target. The scattering process is 
illustrated schematically in Fig. 1(a). The total momen- 
tum k transferred through the scattering of the photon 
is 

k = ki - k2 = q + p„ , (2) 

where the second equality gives the momentum conser- 
vation condition, which is shown also in Fig. l(b)^. 

Care is needed in formulating the energy conservation 
condition. In the standard treatment of the Compton 
scattering process, one assumes independent electrons 
with various one-particle energies. However, the preced- 
ing discussion makes it clear that the general interacting 
system is more naturally characterized via the quantum 
number n of the ionized target. Therefore the relvant 
binding energy i?^"'' is 

eI""^ ^ EniN - 1) - EoiN) (3) 

in terms of the ground state energy of the N-particle sys- 
tem and that of the (A^— 1) particle ionized target. In the 

(n) 

one-particle approximation, Ej^ will correspond to the 
energy of the orbital from which the outgoing electron is 
ejected. The energy conservation then yields 

L^i - W2 = Ej^"^ + . (4) 

In addition to energy and momentum, the total spin is 
also conserved in the scattering process. If the target is 
initially in an 5 = state, then the final state will also 
be a singlet, while for a magnetized sample in a triplet 
state, the final state can be a doublet or a quartet. 

As shown in Refpi^, in the high energy transfer region, 

(n) 

where uji — uj2 » E^ ' , the interaction between the 
electromagnetic field and the target can be approximated 



where A is the vector potential of the field. Using this 
form of the interaction, the expression for the double dif- 
ferential Compton scattering cross section in the lA was 
obtained in Refi^^. The lA corresponds in effect to mod- 
eling the scattering process as an elastic collision between 
a photon and an electron of a particular momentum with 
the target being represented by a distribution of such in- 
dependent electronic states. 

More relevant for our purposes is the treatment of 
Refspi^iifi. These authors obtain the cross section for the 
elementary Compton scattering process (in which the ion 
is left behind in a specific quantum state) from a many- 
body molecular system in the nonrelativistic approx- 
imation of Eq. lOl . It is natural to refer to such a cross 
section as a partial triple differential scattering cross sec- 
tion (PTDSC), since the total triple differential scatter- 
ing cross section (TTDSC) is obtained by summing the 
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PTDSC over final states (see Section HlBl belowl . The 
expression for the PTDSC isiiifi 



cos 



(6) 

where 5 denotes the Kronecker symbol (here and else- 
where in this article) and reflects the energy conservation 
law in Compton scattering, tq = e^/mc^ is the classical 
electron radius, 6 is the scattering angle, and 



M(") *^"^(xi....x^) I "^exp(zkr,) | *,(x 



N 



•Xat) > 



(7) 

is the transition matrix element calculated with N- 
electron wave functions of initial and flnal states of the 
target. Note that expression (jsj for the PTDSC assumes 
an implicit summation over the vibrational states within 
the framework of the Born-Oppenheimer approximation; 
in any event, it will be difficult to resolve vibrational lev- 
els in the Compton scattering regim o^^i^^ . 

The antisymmetry of the many electron wavefunction 
implies that the contribution from each term in expres- 
sion Q is the same. Therefore, we may replace the sum- 
mation over V by the last term, which corresponds to the 
ejection of the A^-th electron, yielding 

]V/(") < ^'^"^(xi....XAr) I exp(ikrAr) | ^',(xi. 



(8) 

We assume now that the initial state possesses a total 
spin 5* = 0, as is the case in most nonmagnetic materials. 
The final state will then be a singlet state (due to spin 
conservation in the Compton scattering process) and the 
associated antisymmetric singlet wave function can be 
represented as 

^'^"^(xi....XAr) = yl--i=[^'„Q(xi....X7V-l)V'p/3(xA') 

-*n/3(xi....Xjv_l)V'pQ(xjv)] • (9) 

Here i/'pcr(xAr) is the wave function of the ejected electron 
with momentum p„ and spin projection a that can accept 
only two values: a = 1/2 denoted by a and tr = — 1/2 by 
/?. ^na and "^nfj are the two doublet components of the 
(TV — l)-electron ionic wave function in the n-th quantum 
state (in the one-electron picture, this describes an ion 
with a hole in the n-th shell). The ionic states are the 
eigenstates of the (TV — l)-electron Hamiltonian. A is an 
antisymmetrization operator given by 



A 



1 



N-l 



(1 - ^ P^n) 



(10) 



where the permutation P^n transposes the ejected TVth 
electron with the vi\i electron in the ion. 

Substituting the final state wave function ^ into Eq. 
and invoking the condition of strong orthogonality 



no. (xi....XAr_i)V'p;3(xAr)- , . 

*„,5(xi....Xjv-l)'0pQ(xAr)] |*i(xi....XAr)) ^ 

)fDr all i, we obtain 



M(") = < [«'„„(xi....xw_i)V'p/3(xw) 

-*„/3(xi....XAr_i)V'pa(xjv)] I exp(zkrjv) 
I *j(xi....xjv) > . 

(12) 

This can be represented in terms of the Dyson spin- 
orbital922i22iiaiiiii^, defined by 

g„(xAr) = \fN j *„(xi....xjv-i)**o(xi....XAr) dxi....(ix 

(13) 

where the integration over dx^ includes a summation over 
the spin coordinates. The Dyson spin-orbitals (7„(xAr) 
may thus be thought of as generalized overlap amplitudes 
between the ground state and the ionized states of the 
many body system. They naturally appear in the spec- 
tral resolution of the one-particle Green functior42iii, 
and have been exploited successfully in some studies of 
ionization of atomic and molecular systems by electro- 
magnetic radiation or fast electrons^^i^Si^. Note that 
in general Dyson orbitals do not form an orthonormal 
set. Some authors^ define Dyson orbitals without the 
pre-factor of ^/N . The Dyson spin-orbital with the spin 
projection a may be written in terms of the spin function 
cr(C) as 

5„(xAf) = 5„(rjv, cr(Cjv)) = 5„(rAr)CT(C^) . (14) 
The wave function of the ejected electron similarly is 



V'p„q(xjv) = V'p„(r7v)a(CAf) 



(15) 



Introducing definition Ijl3|l into Eq. II12II and perform- 
ing spin integration, we obtain a compact general expres- 
sion for the transition matrix element 



M(") = V2j 5„(r)exp(zkr)0;jr) dv 



(16) 



In the region of large energy transfer, Eq. ifTfit pro- 
vides an exact expression for the matrix element in terms 
of the Dyson orbital (7„(xjv) and the wave function of the 
ejected electron in the potential field of the ion. Electron 
correlation effects enter through g„(xAr) and can be in- 
cluded in any particular scheme to the extent to which 
these are incorporated in the computation of this quan- 
tity. In general, Dyson orbitals can be expanded into 
linear combinations of Hartree-Fock or other one-particle 
wave functions. In the so-called diagonal approximation, 
the Dyson orbital is equal to the square root of the pole 
strength times the HF orbitai^i^, and can be calculated 
using special code^ implemented into the Gaussian-98 
program suited; see also Refi^. 
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Under the condition uji — lu2 >> Eb^\ the wave func- 
tion of ejected electron (r) may be approximated as 
a plane wave 



'^P„(r) 



(27r)3/' 



■ exp(ip„r) 



(17) 



allowing the transition matrix element to be expressed 
via the Dyson orbital (?n(q) in momentum space 

^'"^ = 5^ / 9n{Y^ exp(zqr) dr = V2 5„(q), (18) 

Here, q = k — p„ is the momentum transferred to the 
ion. Since the ejected electron is considered as being free 
(with energy p^/2m), the absolute value of the vector p„ 
is completely determined by the energy conservation law 
and is equal to 



select the specific quantum state n involved in the scat- 
tering process through the energy conservation condition 
Q. Another approach, followed more recently by Itoh 
and collaborators, is to fix the position of the electron as 
well as the photon detector (i.e. the angles 9 and 9e in 
Fig. 1(a)), but energy analyze both the scattered photon 
and the recoil electron in coincidenc92SiS£. 



B. Summation over final states 

The TTDSC is obtained from Eq. igll by summing 
over the available final states n as^ 



Pn = \l 2m (Wi - W2 - £^h"'') 



(19) 



The direction of the vector p„ is undetermined so that 
only the maximum and minimum values of the vector q 
are constrained as follows 



\k~Pn\ <q<k+Pn . 



(20) 



In this sense, vector q involves an impHcit dependence 
on the index n. 

Using Eqs. ifTTjl and we obtain the final expression 
for the PTDSC where the ion is created in a definite 
electronic state n as 



dujodVlodD,,- 



= rl{l+cos'e)'^\g,M)\''K^i-^2-Er-^) 



UJl 



2m' 
(21) 



Note that here the independent particle approximation 
(IPM) is not invoked. In the IPM, the Fourier component 
of the Dyson orbital in Eq. Ij21|l reduces to the Fourier 
component of the Hartree-Fock or the Kohn-Sham orbital 
from which the electron is removed. Moreover, aside from 
the use of the plane wave form II17II for the ejected elec- 
tron wave function, expression Ij21|l does not invoke the 
impulse approximation. 

The determination of the PTDSC and 5ri(q) in Eq. 
1I2TI1 requires measurements of the angular and energy 
characteristics of both the scattered photon and the 
ejected electron taken in coincidence. In order to under- 
stand the relevant experimental geometries, it is help- 
ful to refer to the momentum conservation condition de- 
picted in Fig. lb above. Kaplan and Yudini^ suggested a 
scheme in which the characteristics of the outgoing pho- 
ton beam (i.e. the angle 9 and energy lu2 in Fig. 1(a)) 
are fixed, but the angle 9^ of the ejected electron is varied 
to access different q- values. The fixed value of 0^2 should 
be selected near the peak of the Compton \m^^^ (i.e. 
close to the value given by the Compton formula for the 
scattering from a free electron at rest) . By measuring the 
energy of the ejected electron then, one can, in principle. 



d^n 



duJ2d^2d^e 



cos 



UJl 



(") 



Pn , 

2ni' 



(22) 



UJ2» eI"\ the 



In the high-energy transfer region, uii 
binding energy in the Kronecker (5-function on the right 
hand side may be neglected, so that the absolute value 
of momentum p„ in Eq. ifTfljl becomes independent of n 
and the summation over n simplifies to yield 



J2\g4q)\'S{u;i^u;2-Ei 
n 

S{u;i-u;2-^)Y.\gn(.cD\', 



(23) 



where p„ is replaced by Pe to emphasize that the mo- 
mentum of the outgoing electron is independent of state 
n. 

The sum of |(?n(q)P over all occupied states n 
can be expressed via the one-particle reduced density 
matrisSiiS2i5i for an TV-electron system defined as 

ri(r;r') ^Nj «'(xi....xjv-i, r, C)**(xi....XAr_i, r', C) 
X(ixi....c?xjv-idC- 

(24) 

In the momentum space 

ri(q;q) = J ri(r,r')exp(fq(r-r')) dr dv' . 

(25) 

We now recall the following decomposition of the one- 
particle reduced density matris^ 



ri(r;r') = 5]g„(r)5„(r')* . 



(26) 



Substituting this decomposition into Eq. Ij25|l . one ob- 
tains 

^ |5„(q)P = (27r)3ri(q; q) ^ (27r)V(q) • (27) 
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Thus, in the high-energy transfer region, the TTDSC is 
directly related to the 3D EMD as follows 

. (2.)Vg(lW.)^p(q).(..-..-ii) . 

(28) 

It is this TTDSC that is measured in the (7, 67) exper- 
iments by Bell and collaboratorsSSi^SiSi, issues of exper- 
imental resolution notwithstanding. Note however that 
there is an interesting difference in the way the momen- 
tum density factor p occurs in Eq. II28II compared to the 
analytical expressions employed by RefsiS2i22iSi. In our 
case, the EMD (p(q)) is sampled at the momentum 

q = k - Pe , (29) 

which is the momentum transferred to the ion, whereas in 
the cross section of Refsi^^i^i^i^^i^i?'^^?, the EMD involved 
is p(p), where p is the initial momentum of the elec- 
tron before scattering. The reason is that the study of 
j^g£g|22,23^ is based on the formulae of Ribberfors54 for 
double differential cross sections in the lA. As already 
noted, in the lA the scattering is the same as for free elec- 
trons, but weighted with the probability with which the 
plane-wave state of momentum p occurs in the ground 
state. For a system of bound particles, this picture does 
not constitute a useful starting point, and our more gen- 
eral treatment indicates that the quantity that occurs 
naturally is the momentum q transferred to the ion. Nev- 
ertheless, in the lA, the two pictures are equivalent be- 
cause in this regime, from momentum conservation, one 
obtains: 

p -h ki = Pe -t- k2 (30) 

or equivalently 

P = Pe-k. (31) 

Comparing Eq. II29|I and H31|l , we see that in the high en- 
ergy transfer limit, q and p differ only by direction (they 
are opposite), as the outgoing electron loses all memory 
of the bound state it came from. The range of q in for- 
mula Ij2()|l becomes 

< g < 2fc. (32) 

Thus, the maximum momentum transferred to the ion 
is given by 2k. In the lA, 2k may be interpreted as the 



highest momentum of an electron in the initial system 
that can be ejected. 



III. SUMMARY AND CONCLUSIONS 

We start by considering the elementary process in- 
volved in Compton scattering, namely, the scattering 
of an incoming photon from the ground state of an iV- 
electron target to yield a final state containing a singly 
ionized target with (iV— 1) electrons in a specific quantum 
state together with an outgoing photon and an ejected 
electron. The associated PTDSC (partial triple differen- 
tial scattering cross section) is obtained rigorously with- 
out resorting to the approximations inherent in either the 
lA or the IPM. It is shown that the PTDSC can be ex- 
pressed in terms of the Dyson orbitals, which give the 
overlap between the wave function of the ground state 
of the TV-electron initial system with the (N — 1) elec- 
tron ionized final state wave function. The TTDSC (to- 
tal triple differential scattering cross section) is then ob- 
tained by summing over final states, which is equivalent 
to summing over the occupied Dyson orbitals. Interest- 
ingly, in our general treatment, the momentum that plays 
a fundamental role in the formula for the cross section is 
the momentum q transferred to the ion in the scattering 
process and not the momentum p associated with the 
electronic system as is the case in the lA based treatment. 
We show how in the limiting case of the high energy 
transfer regime, our formalism reduces to the standard lA 
description. Although our treatment is non-relativistic, 
extension to the relativistic case is straightforward by us- 
ing relevant results in the literature^^i^.Our formalism 
will permit the analysis and interpretation of electronic 
transitions in correlated electron systems via IXS beyond 
the constraints of the lA and the IPM. 
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Figure Caption 

Figure 1 : (a) Schematic diagram of the elementary scat- 
tering event involved in the Compton scattering process. 
The incoming photon scatters from the target to produce 
an outgoing photon and an electron and leaves the singly 
ionized target (not shown) in a definite quantum state de- 
noted by index n. The notation for kinematic variables 
is obvious, (b) Momentum conservation in the process of 
Fig. 1(a), where q is the momentum transferred to the 
ionized target. 



